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1 Introduction 

Energy storage and attitude control are accomplished with two separate devices on present spacecraft. 
Batteries are typically used to store and supply electrical energy produced by photovoltaic cells; however, 
batteries are quite massive, and battery life often limits the life of a spacecraft. Spinning cylinders, such as 
reaction wheels and control moment gyroscopes, are often employed to control orientation without expending 
propellant. Although these devices possess a great deal of rotational kinetic energy, there are no provisions 
to convert it back into electrical form. Mechanical flywheel systems are an attractive alternative to batteries. 
Their longevity is superior, and they are less massive than batteries; moreover, flywheels can simultaneously 
store energy and control attitude, making it possible to reduce spacecraft mass even further. 

Several advanced technology experiments for the International Space Station (ISS) , including the Attitude 
Control and Energy Storage Experiment (ACESE) led by NASA Lewis Research Center, are currently 
receiving Phase B funding from the Engineering and Research Technology Program at NASA Johnson Space 
Center (JSC). ACESE will advance the development of aerospace flywheels by demonstrating the integration 
of energy storage and attitude control into a single system aboard a working spacecraft. 

The experiment consists of two counter-rotating rotors placed in vacuum housings, and levitated with 
magnetic bearings. Motor-generators will connect the rotors to the existing electrical power system so that 
they can store energy when it is available from the photovoltaic arrays, and supply energy when it is needed. 
Each rotor, made up of a metallic hub and a rim of composite material, will be approximately 11 inches in 
diameter and 13 inches in length, and spin at angular speeds ranging between 15,000 and 50,000 rev/min. 
ACESE is designed to store the same amount of energy as two batteries; it will be attached to the Station’s 
starboard outboard truss at an unoccupied site for a battery charge-discharge unit and two batteries, and 
operate in combination with existing batteries. The primary objective of the experiment is to demonstrate 
energy storage, whereas the secondary objective is to use flywheels to exert torque on the Station, and show 


1 



measurable evidence that torque has been applied. In the interest of simplicity, the time history of the 
magnitude of the torque will be predetermined rather than governed by a feedback control scheme, and the 
existing Station attitude control system will not be modified to account for presence of the flywheels. The 
design of ACESE hardware is simplified considerably by omitting devices for measuring directly the torque 
applied by the experiment to the Station. Measurements can, instead, be obtained indirectly by examining 
telemetry of CMG torque and momentum. The resolution of CMG momentum measurements is expected 
to be approximately 100 ft-lbf-sec, so changes greater than this value can be used to furnish evidence of 
flywheel torque. 

This report presents results of simulations performed to study the behavior of the Station’s attitude, 
and the Control Moment Gyroscopes, when torque is applied by the flywheels in pursuit of the secondary 
objective, or as a result of a malfunction that occurs while storing or discharging energy. Sec. 2 contains a 
description of the simulation program and the model of the flywheel system, a derivation of the algorithm 
for calculating flywheel motor torques, and information regarding mass properties used in the simulations. 
A discussion of CMG momentum capacity and limits placed on momentum produced by Station equipment 
appears in Sec. 3 so that it may serve as a basis for evaluating the results associated with the secondary 
objective, presented in Sec. 4, and the results of simulated malfunctions given in Sec. 5. 

2 Simulation Description 

The numerical results reported herein were obtained with the Space Station Multi Rigid Body Simulation 
(SSMRBS; see Ref. [1]), a computer program that numerically integrates equations of motion governing the 
behavior of a spacecraft modeled as a rigid, multibody, multi-degree-of-freedom system. 

SSMRBS is used at JSC by the Aeroscience and Flight Mechanics Division, support contractors, and the 
Station Guidance, Navigation, and Control (GN&C) prime contractor primarily to evaluate the performance 
of candidate attitude control schemes for Space Station, and to validate GN&C systems for each stage of 
the Station assembly sequence. Analysts are able to study the interplay of attitude control with other Space 
Station activities such as collection of solar power, radiation of excess heat, remote manipulator operations, 
centrifuge operations, and microgravity experiments. Nominal operations and contingencies are studied, 
including control of a Station with or without an Orbiter attached, Orbiter docking or separation, and 
Station reboost. 

Simulations typically involve, simultaneously, the control of core body attitude with a control moment 
gyroscope momentum manager or a reaction control system (or both) ; independent feedback control of each 
solar array and radiator joint angle; and transportation of payloads with a mobile transporter or remote 
manipulator. Analysis involving the ACESE flywheels has been performed at JSC and, more recently, at 
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the NASA Langley Research Center. 

2.1 Model of the Flywheel System 

All bodies are treated as rigid in the simulations to be discussed; a flywheel rotor is considered to be a right 
circular cylinder whose mass is distributed uniformly. The model of the physical connection between a rotor 
and the outboard truss that is simplest and most appropriate for the study of the Station’s attitude behavior 
is a revolute, or hinge joint whose axis is coincident with the rotor’s longitudinal central principal axis of 
inertia, or spin axis. In such a model of the joint, the spin axis of the rotor is fixed in the outboard truss 
and in the rotor; that is, the model does not account for any changes in the spin axis direction relative to 
the outboard truss that are in fact allowed by the magnetic bearings, nor does it permit any coning motions 
of the rotor. Counter rotation of the two rotors requires that their spin axes be parallel. 

Relative motion between a body and its inboard body can be brought about in the simulations by means 
of a motor; the masses of the two bodies are considered to include the mass of the motor parts. The inboard 
body exerts torque on the outboard body through the motor, and in accordance with the law of action 
and reaction, torque of equal magnitude and opposite direction is applied to the inboard body. An ACESE 
motor-generator is represented by such a motor, whose detailed electrical behavior need not be modeled. 
The focus of the remainder of this section is the calculation of the torque that must be exerted by each 
of two motors in order to apply a specified resultant torque to the outboard truss, and change the kinetic 
energy of the rotor pair at a specified rate. 

We begin by denoting the torque exerted by the outboard truss C on rotor A as the vector T C ^ A , and 
the torque exerted by C on rotor B as T c ! B . The resultant of the two is important in a discussion of the 
Station’s attitude motion, and is given by 

T °/F t T c / a + T c/b (1) 

where F is the system of flywheels formed by rotors A and B. 

The rate of change of rotational kinetic energy, or power of the rotors is an important measure of energy 
storage; thus it is to be specified. The rotational kinetic energy C K F of F in C is given by 

°K f = (°u A ■ \ A ' A ‘ ■ c o J A + c u> B . \ B ' B * • c u, B ) (2) 
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where c uj j4 and c u> B are, respectively, the angular velocities of A and B in (or relative to) C, and where 
and 1 B ' B * are the inertia dyadics of A and B relative to their respective mass centers, A * and B*. 
The power c P F of F in C is defined as the derivative of c K F with respect to time, therefore 


— c k f 
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where all vectors and dyadics are differentiated with respect to time in the same reference frame C, as 
indicated with the notation c d/dt. One may recognize, first, that I' 4/< ' 4 and 1 B / B do not change with time 
in C because the joints connecting A and B to C are parallel to the axes of symmetry of the rotors, and 
second, that the angular accelerations c a A and c a. B of A and B in C are defined respectively by 


c c 

C a A A ±C U A C a B A ± CuJ B 

dt ’ dt 


Hence, c P F can be expressed as 


'P F = (r 4 / A ‘ . CV 4 ) . C U A + ^B/B * 


C B\ C, , B 


By appealing to the angular momentum principle, and assuming that the moment of all forces exerted 
on a rotor about the rotor’s mass center is equal to the torque exerted by C on the rotor, we may write 


T°/A = JA/A* .C a A, C U A x j A/A* . C 4 


t C/B = J B/B* . C a B , C 


UJ B xl B/B*,C U B (6) 


Now, the revolute joint ensures that the angular velocity c uj a is always parallel to a central principal axis 
of inertia of A; therefore, the term c ui A x I' 4 /' 4 • c u i A vanishes. Likewise, the counterpart of this term for 
rotor B also vanishes. We thus obtain 


t c/A = jA/a* m c a A 


T C / B = J B/B* . C a B 


Substitution from Eqs. (7) and (1) into (5) yields 

Cp F _ rpC/A _ C y /1 _|_ rjtC/B _ C^B 

_ ^ T C/F _ T C/B ^ A + T C/B . C^B 

^ T c/F . C M A + t C/B . ^C^B _ C^A'J ( 8 ) 

In order to proceed to express the right hand member in terms of scalars, it is convenient to introduce a 
unit vector A parallel to the rotor spin axes, fixed in rotors A and B, and in the outboard truss C. We then 
define the following three scalars 

rjnCjA rj-iO/A ^ rpC j B ^ rji CjB ^ rj-\C j F ^ rj^O j F ^ (9) 
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and write the angular velocities of A and B in C as 


A A c 
L U! ' = M4A, 


C B A t 

u> = u B A 


(10) 


where u A and u B are angular speeds in C of A and B , respectively. After substituting from Eqns. (9) 
and (10) into (8), forming the indicated scalar products, and solving the resulting expression for T C ^ B , one 
obtains 


rpC/B 


CpF _ T C, IP u A 

u B -U A 


( 11 ) 


In addition, the scalar product of Eq. (1) with A yields, after rearrangement, 


rpC/A __ rpC/F _ rpC/B 


( 12 ) 


With the c P F and T C ' F specified, and the rotor angular speeds u A and u B known, Eqs. (11) and (12) are 
used 1 , in order, to calculate the A-measure numbers of the torques to be applied by motors to B and A. 
With °P F > 0, the flywheels are charging; that is, rotational kinetic energy in C of F is increasing. Torque 
is applied by F to C in the direction of A when T c ! F < 0. 

When ACESE is to be used solely to store and supply energy, without exerting torque on the Station, 
T c ! f is set to zero. In addition, the angular speeds of the rotors should begin with identical values and 
opposite signs; that is, u A = —u B . In view of Eq. (11), T c ! B = c P F /(2u B ), and, in view of Eq. (12), 
T c ' a = —T c / b . Under these conditions, the absolute values of u A and u B will remain identical. 


2.2 Mass Distribution of Station Bodies, Rotors 

The Space Station is modeled as a collection of 15 rigid bodies fastened together: a core body, outboard truss 
structures, solar arrays, radiators, and so forth. The mass distribution for each body, and the topography of 
the spacecraft are as described for configuration c080_13a in Ref. [2], Two additional bodies attached to the 
starboard outboard truss represent the flywheel rotors. The masses of A and B are given by m A = rn B = 2.17 
slugs, and the central inertia dyadics are given by T 4 /" 4 = \ B / B = 0.294 CiCj + 0.222 AA + 0.294 6363 
slug- ft", where c,\ and £3 are any two unit vectors perpendicular to each other and to A. 


3 CMG Momentum Capacity, Payload Momentum Limits 

Evaluation of the following simulation results is facilitated by knowing the angular momentum capacity of 
the Station’s Control Moment Gyroscopes (CMGs), and limits placed upon angular momentum that can be 
produced by payloads or experiments during the time that microgravity experiments are being conducted. 

1 Eqns. (11) and (12) are in agreement with Eqs. (7) and (8) in notes provided by Michael Oshima at Boeing (Seattle) 
on Sept. 2, 1997. 
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The magnitude of central angular momentum in an inertial reference frame that each CMG can possess 
is limited to 3500 ft-lbf-sec; when all four CMGs are operating, the magnitude of the resultant of the angular 
momenta can reach 14000 ft-lbf-sec. However, the Station’s reaction control system jets will expend pro- 
pellant (and disturb any microgravity experiments in progress) to desaturate the CMGs if their momentum 
reaches 13000 ft-lbf-sec. In comparison, one flywheel rotor possesses approximately 1160 ft-lbf-sec of angular 
momentum at high speed (0.222 slug-ft 2 x 50000 rev/min); that is, about A of the momentum of a CMG. It 
is useful to keep this in mind when reading Secs. 4 and 5. 

During nominal operations aboard the Assembly Complete configuration, equipment on the United States 
On-orbit Segment (which includes the outboard truss structures) is not permitted to impart more than 1600 
ft-lbf-sec (see Ref. [3]) to the Station while microgravity experiments are being conducted. Although this 
limit is not strictly applicable to the 13a configuration, or to contingencies, it nevertheless furnishes a useful 
measure for assessing the results reported in Sec. 5. 

4 Flywheel Torque Assist for CMGs 

The goal of using ACESE to exert torque on the Station, and demonstrate measurable evidence thereof, gives 
rise to questions of the direction and magnitude of the torque that will be best for this purpose. In terms of 
the symbols already introduced, we wish to know how should A be oriented in C, and what mathematical 
form should T C ! F take. 

As it turns out, the orientation of the rotor spin axes in C is dictated more by the geometry of the 
ACESE enclosure and the site to which it is to be attached, than by considerations of the secondary goal. 
However, it is fortuitous that the geometrical constraints lead to a direction of A that facilitates attitude 
control. The rotors are packaged most easily with A in the same direction as 62 , a unit vector fixed in the 
outboard truss and parallel to the axis of the revolute joint (known as the alpha joint) that connects the 
outboard truss to the Station core body. Unit vector c 2 is thus fixed also in the core body, in what is called 
the Station “y” direction, leading to analytical and operational simplifications regarding the direction of the 
torque applied by rotors A and B to the Station. 

A sinusoidal form for T C ! F is suggested by two aspects of the flywheel experiment. First, as the Station 
passes from sunlight into the Earth’s shadow, and back into sunlight, the flywheels will alternately charge with 
electricity from the solar arrays, and discharge their energy into the power system. Second, exerting torque 
with the flywheel pair will require their spin speeds to become unequal, and this process must eventually be 
reversed so that the flywheel pair will again be in a position to store energy without exerting torque on the 
Station. 

If possible, torque should be applied by the flywheels in such a way as to assist the CMGs in their job of 
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controlling the orientation of the Station core body, keeping it in “torque equilibrium attitude” . Although it 
may be possible to obtain measurable evidence of applied torque with an arbitrary time history of T c / F , it 
makes more sense to anticipate the torque to be exerted by the CMGs on the core body, and attempt to supply 
part of this with the flywheels. Simulation results for Station configuration 13A, with T C ! A = T C ^ B = 0, 
are shown in Fig. 1 for a time interval during which attitude motion has reached steady-state; the torque 
exerted by the CMGs on the core body is projected into unit vectors bi, b 2 , and b 3 fixed in the core body 
and parallel to the Station “x” , “y” , and “z” directions respectively. The alpha joint that attaches the core 
body to the outboard truss keeps b 2 equal to c 2 (and therefore to A) at all times. Inspection of Fig. 1 reveals 
that the b 2 measure number of CMG torque is nearly a simple sinusoid, whereas the other two curves can 
not be represented well with this simple mathematical form. The orbital period of the Station is in this case 
approximately 5560 sec, so the sinusoid has a frequency of two cycles per orbit. 


c080_1 3a DAC7; with ACESE flywheels hS - li it Nm 

spin axis parallel to outboard truss y 
no flywheel torque exerted on core body s 



20.0 21.0 22.0 23.0 24.0 25.0 25.0 27.0 28.0 29.0 30.0 W 

Time (Sec) *10 3 < 

CD 



Figure 1: CMG Torque, No Flywheel Motion 


Consequently, the form chosen for T C ! F is 

T c/F = Tsm2u;(t-t*) (13) 

where u> is the frequency at which the Station orbits the Earth, and t* is a particular value of the time 
t at which the b 2 measure number of CMG torque is zero, and has a positive slope. The angular speeds 
ua and ub of the rotors should remain, approximately, less than 50,000 rev/min, and greater than 15,000 
rev/min; this requirement, together with the value of the axial central principal moment of inertia of the 
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rotors, effectively limits the absolute value of r to 0.75 ft-lbf (1.017 N-m). The second plot from Fig. 1 is 
redisplayed in Fig. 2 with a solid curve, and a plot of T F ! C = —T C ! F is displayed with a dashed curve, 
where T C ! F is given by Eq. (13) with r = —0.75 ft-lbf, u; = 0.0011295 rad/sec, and t* = 20860 sec. Prior 
to t = t*, t = 0. 


CMG Torque (no f/w motion), Flywheel Torque (ft-lbf) 



Figure 2: CMG Torque, Flywheel Torque (T F / C ) 


A simulation of using the flywheels for simultaneous assistance of the CMGs and energy storage is 
performed, where T C ! B and T C ^ A are calculated according to Eqs. (11) and (12) for t > t*. c p F is 
specified as 2200 watts during the sunlit portion of an orbit, and as —1.57 X 2200 watts during eclipse. (In 
these particular simulations, the ratio of time spent by the spacecraft in sunlight to that spent in eclipse 
is 1.57.) T f /° is specified as indicated in Fig. 2. The results are exhibited with a solid curve in the 
upper plots of Figs. (3)— (5) , and compared to results of the simulation in which flywheel torque is absent 
(T C / A = t c / b = o), shown with a broken curve. The lower plot shows the difference between the two 
curves in the upper plot; specifically, the values displayed with the solid curve (flywheel torque present) are 
subtracted from the values displayed with the broken curve (flywheel torque absent) . 

The reduction in the ba measure number of torque exerted by the CMGs as a consequence of assistance 
from the flywheels is illustrated in Fig. (3) over the interval 30000 sec < t < 40000 sec. The decrease, 
quantified in the lower plot, is seen to have the form of Eq. (13), with an amplitude of very nearly 0.75 ft-lbf. 

The reduction in CMG torque is reflected in a smaller magnitude of the resultant of central angular 
momenta in inertial space of the four CMGs, as shown in Fig. 4. The difference is displayed in the lower 
plot; a positive value in this case indicates a reduction in momentum magnitude, which periodically reaches 
200 ft-lbf-sec. 

The CMG momentum management scheme used in these simulations places greater emphasis on control- 
ling orientation than on minimizing CMG momentum; therefore, one expects to see only small differences 
in the roll, pitch, and yaw of the core body with respect to, say, a local-vertical-local-horizontal reference 
frame. Such expectations are in fact borne out; of the three orientation angles, the largest difference occurs 
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in pitch and is less than 0.18 deg, and is transient in nature as shown in Fig. 5. 

The foregoing differences in CMG torque and momentum should be in evidence in telemetry transmitted 
from the Station, providing measurable confirmation of the sinusoidal torque exerted by the flywheels. 
Differences in core body attitude are less likely to furnish a clear demonstration of the experiment’s success. 




nowheels — wheels (deg) *10 -z Ivp *10° nowheels — wheels (ft— Ibf— sec^io 1 







5 Energy Storage Contingencies 

Ideally, the storage and discharge of energy takes place as the rotors are counter rotating ( c u> A = — 1 a> B ), 
and no torque is applied by the flywheel pair to the outboard truss C . However, malfunctions can occur 
in which the rotors no longer counter rotate, and T F ! C differs from 0. One such situation involves having 
to bring a rotor to an abrupt stop while the other rotor continues to spin. A second condition of this kind 
occurs when energy storage must be performed with one rotor because it is not possible to operate the second 
rotor. Both of these situations are examined in what follows. 

5.1 Emergency Shutdown of One Rotor 

An abrupt stop of one rotor is simulated by starting with initial values ug(t = 0) = —u^(t = 0) = 50000 
rev/min, setting T C ^ A = T C ^ B = 0 for 0 < t < 20860 sec, and then making T C ! B = —2 ft-lbf for 
20860 < t < 21441 sec. In this way B is brought to rest in C, while the angular speed ua remains equal 
to —50000 rev/min. As in the previous section, results are compared to those of the simulation in which 
T c / a = T°! b = 0 at all times, and UA(t = 0) = us(t = 0) = 0. Several parameters are examined, 
and the greatest differences observed in each are reported in Table 1. “Attitude, wrt LVLH” refers to the 
orientation of the Station core body relative to a local-vertical-local-horizontal reference frame, described 
with a body-three, 3-2-1 (yaw, pitch, roll) rotation sequence. The b]-b 2 -b >3 measure numbers of the resultant 
of the four CMG central angular momenta, and the resultant CMG torque, are examined. In addition, the 
magnitude of the CMG momentum is shown in the final column. Time histories of two parameters are shown: 
Fig. 6 displays the core body pitch angle, and Fig. 7 shows the b 2 -measure number of CMG momentum. 
The differences in attitude and CMG torque are not significant, and the difference in CMG momentum 
magnitude is within the limits set for payload angular momentum during microgravity operations aboard 
Assembly Complete (see Sec. 3). 


Table 1 : Effects of Rotor Shutdown 



yaw 

pitch 

roll 


Attitude, wrt LVLH (deg) 

0.023 

0.350 

0.035 


• bf 

•b 2 

' b 3 

magnitude 

CMG Angular Momentum (ft-lbf-sec) 

75 

1600 

100 

850 

CMG Torque (ft-lbf) 

0.12 

3.3 

0.19 

- 
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c080_ 1 3a DAC7; with ACESE flywheels IK Elh 9 .lt, 

spin axis parallel to outboard truss y 
overspeed study: —2 ft— Ibf on wheel B 



Figure 6: Core Body Pitch Relative to LVLH 



Figure 7: ba measure number of CMG momentum 
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5.2 Energy Storage Using a Single Rotor 


Should one rotor become incapacitated, it may be possible to demonstrate energy storage at a reduced level 
with the remaining rotor. Supposing that A is regarded as the malfunctioning rotor, the event is simulated 
by setting T C ! F to 0 in Eq. (11), and by specifying c P F as 1100 watts when the spacecraft receives sunlight, 
and —1.57 X 1100 watts when the spacecraft is in darkness. Instead of calculating T c ! A according to Eq. (12), 
it is defined to be 0 throughout the simulation. 

The particular initial values of orbit parameters and position of the sun in these simulations make it 
convenient to choose us{t = 0) = 38000 rev/min; thereafter, ub varies between 51000 and 19000 rev/min. 
Since A is to be at rest in C, ua is chosen as 0 at t = 0, and remains exceedingly small for the duration of 
the simulation. 

As before, results are compared to those of the simulation in which p c l A = pC/B = q at all times, 
and the greatest differences observed are reported in Table 2. Time histories of core body pitch and the 
b 2 -measure number of CMG momentum are shown in Figs. 8 and 9 respectively. The differences seen in 
attitude and CMG torque are not significant, and the difference in CMG momentum magnitude is within 
the limits set for payload angular momentum during microgravity operations aboard Assembly Complete 
(see Sec. 3). 


Table 2: Effects of Single- Rotor Energy Storage 



yaw 

pitch 

roll 


Attitude, wrt LVLH (deg) 

0.010 

0.140 

0.019 


•bj 

• b 2 

• t> 3 

magnitude 

CMG Angular Momentum (ft-lbf-sec) 

55 

650 

100 

300 

CMG Torque (ft-lbf) 

0.12 

1.2 

0.12 

- 
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single wheel operation in power mode, wheel B 
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Figure 8: Core Body Pitch Relative to LVLH 
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Figure 9: t >2 measure number of CMG momentum 










6 Conclusions 


The secondary objective of the ACESE experiment, to contribute torque for Station attitude control while 
storing or discharging energy, can be achieved by matching the phase of the sinusoidal torque exerted by 
CMGs when attitude motion reaches steady state. Although the amplitude of the torque exerted by the 
flywheels is much less than that of the CMGs, a reduction in CMG torque and momentum should be evident 
in telemetry of those parameters. 

Two ACESE contingencies were studied: the abrupt stop of one rotor while the other rotor continues 
to spin at high speed, and energy storage with one rotor instead of a counter rotating pair. In each case, 
Station attitude and CMG torque were not affected significantly, and CMG momentum remains within limits 
established for microgravity operations aboard the Assembly Complete station configuration. 
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